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The quantum capacity of a memoryless channel is often used as a single figure of merit to 
characterize its ability to transmit quantum information coherently. The capacity determines 
the maximal rate at which we can code reliably over asymptotically many uses of the channel. 
We argue that this asymptotic treatment is insufficient to the point of being irrelevant in the 
quantum setting where decoherence severely limits our ability to manipulate large quantum 
systems in the encoder and decoder. For all practical purposes we should instead focus on 
the trade-off between three parameters: the rate of the code, the number of coherent uses 
of the channel, and the fidelity of the transmission. The aim is then to specify the region 
determined by allowed combinations of these parameters. 

Towards this goal, we find approximate and exact characterizations of the region of al¬ 
lowed triplets for the qubit dephasing channel and for the erasure channel with classical 
post-processing assistance. In each case the region is parametrized by a second channel 
parameter, the quantum channel dispersion. In the process we also develop several general 
inner (achievable) and outer (converse) bounds on the coding region that are valid for all 
finite-dimensional quantum channels and can be computed efficiently. Applied to the depo¬ 
larizing channel, this allows us to determine a lower bound on the number of coherent uses 
of the channel necessary to witness super-additivity of the coherent information. 


I. INTRODUCTION 

One of the quintessential topics in quantum information theory is the study of reliable quantum 
information transmission over a noisy quantum channel. Here the word “channel” simply refers 
to a description of a physical evolution (by means of a completely positive trace-preserving map 
on density operators). Traditionally one considers point-to-point communication settings where 
a memoryless channel can be used many times in sequence. The sender (often called Alice) first 
encodes a quantum state into a sequence of registers and then sends them one by one trough 
the channel to the receiver (often called Bob). Bob collects these registers and then attempts to 
decode the quantum state. Alternatively, consider a collection of physical qubits that are exposed 
to independent noise. The goal is then to encode quantum information (logical qubits) into this 
system so that the quantum information can be decoded with high fidelity at a later stage. One 
of the primary goals of information theory is to find fundamental limits imposed on any coding 
scheme that tries to accomplish such tasks. 

Following a tradition going back to Shannon’s groundbreaking work [32], this problem is usually 
studied asymptotically: the quantum capacity of a channel is defined as the optimal rate (in qubits 
per use of the channel) at which we can transmit quantum information with vanishing error as 
the number of channel uses goes to infinity. In the context of information storage, the rate simply 
corresponds to the ratio of logical to physical qubits, and the number of channel uses corresponds 
to the number of physical qubits. The quantum capacity of arbitrary channels has been determined 
in a series of works, an upper bound to the capacity shown in [2, 3, 31] and achievability of that 
bound shown in [9, 22, 33]. 

However, in any application of the theory resources are finite and the number of channel uses is 
necessarily limited. More importantly, at least for the near future it appears unrealistic to expect 
that encoding and decoding circuits can coherently manipulate large numbers of qubits. Restricting 
the size of the quantum devices used for encoding the channel inputs and decoding its outputs is 
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FIG. 1: Coding Scheme for entanglement transmission over n uses of a channel M = Ma^b- The systems 
M, M' and M” are isomorphic. The encoder £ = encodes the part M' of the maximally entangled 

state 4>mm' into the channel input systems. Later, the decoder V = Vb^^m" recovers the state from the 
channel output systems. The performance of the code is measured using the fidelity , Pmm")- 


tantamount to considering communication with only a fixed number of channel uses. This then 
raises the question whether an asymptotic approach where this number goes to infinity — which 
has proven to be very successful for the analysis of classical communication systems — is equally 
suitable for the quantum setting. Clearly, what we really want to understand is how well we can 
transmit quantum information in a realistic setting where the number of channel uses and the size 
of quantum devices is limited. The quantum capacity is at most a proxy for the answer to this 
question, and in this article we argue that it is often not a very good one. 

The study of such non-asymptotic scenarios has recently garnered significant attention in clas¬ 
sical information theory [16, 27, 38] as well as in quantum information theory [8, 21, 40, 41]. Here 
we extend these considerations to the setting of quantum communication. 

Outline: The remainder of the paper is structured as follows. In Section II we discuss our main 
results detail. In Section III we formally introduce relevant notation and information measures. In 
Section IV A we derive our converse (outer) bounds and in Section IV B we derive our achievability 
(inner) bounds. Finally, Section V discusses the specific examples presented below as Results 1, 2, 
and 3. 


II. DISCUSSION OF RESULTS 

In this work we focus on codes enabling a state entangled with a reference system to be reliably 
transmitted through the channel. This is a strong requirement: reliable entanglement transmission 
implies reliable transmission, on average, of all non-entangled input states. The coding scheme is 
depicted in Figure 1. We are given a quantum channel J\f = Ma^b and denote by AA®" the n-fold 
parallel repetition of this channel. An entanglement transmission code for AA®"' is given by a triplet 
{\M\, £,!£>}, where |M| is the local dimension of a maximally entangled state is to be 

transmitted over A/’®'^. The quantum channels £ = and V = are encoding and 

decoding operations, respectively. (A more formal treatment will follow in Section IIIB.) With 
this in hand, we now say that a triplet {R, n, e} is achievable on the channel J\f if there exists an 
entanglement transmission code satisfying 

ilog|M|>R and (^ o AA®” o £:)((/)^j\^/)^ > 1 - e . 


( 1 ) 
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Here, R is the rate of the code, n is the number of channel uses, and e is the tolerated error 
measured in terms of the fidelity F. 

The non-asymptotic achievable region of a quantum channel M is then given by the union of 
all achievable triplets {R,n,£}. The goal of (non-asymptotic) information theory is to find tight 
bounds on this achievable region, in particular to determine if certain triplets are outside the 
achievable region and thus forbidden. For this purpose, we define its boundary 


Rjsf{n; e) := max : (R, n, e) is achievable on AA| , (2) 

and investigate it as a function of n for a fixed value of We will often drop the subscript M if 
it is clear which channel is considered. 

To begin, let us rephrase the seminal capacity results [2, 3, 11, 22, 31, 33] in this language. The 
quantum capacity is defined as the asymptotic limit of Rjv{n] e) when n (first) goes to infinity and 
e vanishes. The capacity can be expressed in terms of a regularized coherent information; 


QiM) := lim lim Rj^{n;e) = lim 

e—^0 n^oo i^oo 


/c(AA®^) 


( 3 ) 


where the coherent information R will be defined in Section IV B. This result is highly unsatisfac¬ 
tory, not least because the regularization makes its computation intractable.^ Worse, the statement 
is not as strong as we would like it to be because it does not give any indication of the fundamental 
limits for finite £ or finite n. 

For example, even sticking to the asymptotic limit for now, we might be willing to admit a 
small but nonzero error in our recovery. Formally, instead of requiring that the error vanishes 
asymptotically, we only require that it does not exceed a certain threshold, e. Can we then achieve 
a higher asymptotic rate in the above sense? Surprisingly, the answer to this question is not known 
in general. Recent work [42] at least settles the question in the negative for a class of generalized 
dephasing channels and in particular for the qubit dephasing channel 


{l-'y)p + -fZpZ , (4) 

where 7 G [0,1] is a parameter and Z is the Pauli Z operator. Dephasing channels are particularly 
interesting examples because dephasing noise is dominant in many physical implementations of 
qubits. The results of [42] thus allow us to fully characterize the achievable region in the limit 
n —)• 00 for such channels, and in particular ensure that 

lim (n;e) = Ic{Z^), (5) 

n^oo ' 

independent of the value of e G (0,1). Note also that the regularization is not required here since 
these channels are degradable [ 10 ]. 

Here we go beyond studying the problem in the asymptotic limit and develop characterizations 
of the achievable region for finite values of n. We find inner (achievability) and outer (converse) 
bounds on the boundary of the achievable region. These do not agree for general channels (which 
is unsurprising given the fact that such an agreement has not even been established asymptotically 
for nonzero error), but they do coincide for certain important examples. 


^ An alternative approach would be to investigate the boundary ij\f{n-,R) := max{e : {R,n,e) is achievable}. This 
leads to the study of error exponents (and the reliability function) as well as strong converse exponents. We will 
not discuss this here since such an analysis usually does not yield sufficiently tight bounds for small values of n. 
It is not clear if the limit ^ 00 is necessary for any fixed channel, but it was recently shown that there does not 

exist a universal constant R such that C{Af) < for all channels N [7]. 
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(a) Boundary of the achievable region for different values of n 
(second order approximation). 
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(b) Boundary of the achievable region for e = 5% 
(third order approximation in (6)). 


(c) Comparison of strict bounds with third order 
approximation for e = 5%. 


FIG. 2: Approximation of the non-asymptotic achievable rate region of a qubit dephasing channel with 
7 = 0.1 (see Result 1). 


A. Qubit Dephasing Channel 


The first example is the qubit dephasing channel. Building on recent work that established 
the strong converse for this channel [42], we will show that its non-asymptotic achievable region is 
equivalent to the corresponding region of a (classical) binary symmetric channel. This allows us 
to employ results from classical information theory and establish the following characterization of 
the achievable region for the qubit dephasing channel. 

Result 1. For the qubit dephasing channel with 7 G [0,1], the boundary R(n-,e) satisfies 


R{n; e) 


1 - /i(7) -b 


bh $-!(,)+ !^ + o 


n 


2n 



( 6 ) 


where is the cumulative standard Gaussian distribution, h{'y) := —7 log 7 — (1 — 7 )log(l — 7 ) 
denotes the binary entropy and v(j) the corresponding variance, v{'y) := 7 (log 7 + ^( 7 ))^ + (1 — 

7)(log(l-7) + /i(7))2. 

The expression without the remainder term 0 ( 7 ) is called the third order approximation of the 
(boundary of the) non-asymptotic achievable region. It is visualized in Figures 2 for an example 
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(a) Boundary of the achievable region. 
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(b) Comparison of exact bonnds with third order 
approximation. 


FIG. 3: Approximation of the non-asymptotic achievable rate region (with classical post-processing assis¬ 
tance) of a qubit erasure channel with /3 = 0.25 and error parameter e = 1% (see Result 2). 


channel with 7 = 0.1. In Figure 2(a) we plot the smallest achievable error e as a function of the rate 
R. Here we use the second order expansion without the term 2 ^ log n since it can conveniently be 
solved for e. In the limit n —)> 00 we see an instantaneous transition of e from 0 to 1, the signature 
of a strong converse: coding below the capacity Q{Z^) = 1 — h{'y) is possible with perfect fidelity 
whereas coding above the capacity will necessarily result in a vanishing fidelity. 

In Figure 2(b) we plot the third order approximation in ( 6 ) for the highest achievable rate, 
R{n; e), as a function of n for a fixed fidelity of 95% (i.e. we set e = 5%). For example, this allows 
us to calculate how many times we need to use the channel in order to approximately achieve 
the quantum capacity. The third order approximation shows that we need approximately 850 
channel uses to achieve 90% of the quantum capacity. Note that a coding scheme achieving this 
would probably require us to coherently manipulate 850 qubits in the decoder, which appears to 
be a quite challenging task. This example shows that the capacity does not suffice to characterize 
the ability of a quantum channel to transmit information, and further motivates the study of the 
achievable region for finite n. 

Finally, we remark that the third order approximation is quite strong even for small n. To prove 
this we compare it to concrete upper and lower bounds on R{n;e) in Figure 2(c) and see that the 
remainder term O(^) becomes negligible for fairly small n ~ 100 for the present values of 7 and e. 


B. Qubit Erasure Channel 

Another channel we can analyze in this manner is the qubit erasure channel, given by the map 

£0 : {1-(3)p +P\e){e\ , (7) 

where j3 G [0,1] is a parameter and |e)(e| is a pure state orthogonal to p. Here we investigate 
coding schemes that allow classical post-processing (cpp) between the sender and receiver after the 
quantum transmission (see also Figure 5 in Section HI). We denote the corresponding boundary 
of the achievable region by i?^PP(n;e). Since this includes all codes that do not take advantage of 
cpp, we clearly have R{n]e) < R^^'^{n;e) for all channels. 

Here we can determine the boundary i2^PP(n;e) exactly, again by generalizing [42] and relating 
the problem to that of the classical erasure channel. 
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(b) Exact outer bound for a = 0.0825 and e = 5.5%. 


FIG. 4: Approximate inner and outer bounds on the non-asymptotic achievable rate region for the depolar¬ 
izing channel (see Results 3 and 5). 


Result 2. For the qubit erasure channel with j3 G [0,1], the boundary R‘^PP(n;e) satisfies 

e= ^ . (8) 

Z=n-A:-|-1 ^ 

Moreover, for large n, we have the expansion 

R"PP(n;e) = 1 - fi + t/ ~ + of-) • (9) 

V n \n/ 

The latter expression is a third order approximation of the achievable region, but this time 
the term proportional to ^logn vanishes. In Figure 3 we show this approximation for a qubit 
erasure channel with fi = 0.25 and fidelity 99%. In Figure 3(a) we see that the non-asymptotic 
achievable region reaches 90% of the channel capacity for n ~ 180. Again, this confirms that the 
non-asymptotic treatment is crucial in the quantum setting. In Figure 3(b) we compare the third 
order approximation with the exact boundary of the achievable region in (8). We see that the 
approximation is already very precise (and the term O(^) thus negligible) for fairly small n ~ 50. 


C. Depolarizing Channel 

Another basic channel of interest is the qubit depolarizing channel. It is given by the map 

V^-. {l-a)p+^{XpX+ YpY + ZpZ) , (10) 

O 

where a G [0,1] is a parameter and X, Y, Z are the Pauli operators. For this channel not even 
a closed formula for the quantum capacity Q{'Da) is known, and the non-regularized coherent 
information 


= 1 - h{a) - a log 3 (11) 

is only a strict lower bound on it [12] (where h{a) again denotes the binary entropy). However, 
various upper bounds on the quantum capacity have been established as well [29, 34, 35, 37]. For 
example, in [34] it is shown that Q{'Da) < Q{Za) = 1 — /i(a), the quantum capacity of the qubit 
dephasing channel with dephasing parameter a. Here we extend this result to the non-asymptotic 
setting and find the following outer (converse) bound for the achievable rate region. 
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Result 3. For the qubit depolarizing channel T)a with a G [0,1], the boundary satisfies 

Rvcin;£) < RzM;e), ( 12 ) 

where Rz^ (^5 denotes the boundary of the achievable rate region for the qubit dephasing channel 
with dephasing parameter a as in Result 1. 

Clearly this allows us to recycle the bounds in Result 1 and use them as outer bounds for the 
achievable region. This is done in Figure 4 for two example channels. In Figure 4(a) we plot the 
second order approximation of the outer bound for a depolarizing channel with a = 0.05 and 99% 
fidelity. We can see that in order to exceed the coherent information, we will need to code for at 
least Nq = 738 channel uses. This indicates that the question of whether the coherent information 
is a good or bad lower bound on the asymptotic quantum capacity is not practically relevant as 
long as we do not have a quantum computer that is able to perform a decoding operation on many 
hundreds of qubits. We also show a second order approximation for a general inner bound which 
is given in Result 5 below. 

In Figure 4(b) we examine a channel with parameters a = 0.0825 and e = 5.5%. Instead of 
using an approximation for the outer bound we use the exact outer bound to give the answer 
(it is 42) to the question of how many channel uses we need at minimum to exceed the coherent 
information. However, note that this does not give us any indication of what code (in particular if 
it is assisted or not) one would need to use for this purpose. 


D. General Bounds 

We have so far focused our attention on two specific (albeit very important) examples of chan¬ 
nels. However, many of the results derived in this paper also hold more generally. For example, we 
find the following outer (converse) bound for coding schemes that allow classical post-processing. 

Result 4. For any quantum channel M, the boundary R'^PP(n;e) satisfies 

R^PP(n;e)<-log/(AA®%e), (13) 

where f{M,e) is the solution to a semidefinite optimization program. Moreover, if M is covariant 
we find the asymptotic expansion 

R‘=PP(n;e) < RZterin-,^), with RZZer{n-,£) = , (14) 

V n V ^ / 

where the Rains information, and its variance, are defined in Theorem 3. 

The semidefinite optimization program f{J\f,e) (see Section IV A for details) is similar in spirit 
to the metaconverse for classical channel coding by Polyanskiy et al. [27], formulated as a linear 
program by Matthews [23].^ Note that the bound (14) is tight up to the second order asymptotically 
for the qubit dephasing channel (Result 1) and the erasure channel with classical post-processing 
assistance (Result 2). However, in the generic covariant case the bound is in general not expected to 
be tight even in first order. Moreover, if the channel is not covariant we cannot find any non-trivial 
outer bounds on the non-asymptotic achievable region that allows for an asymptotic expansion in 
the above sense. 

Moreover, an inner (achievability) bound of the form shown in Result 1 also holds generally for 
all quantum channels. 


® For quantum coding, Matthews and Leung [20] also give semidefinite optimization program lower bounds on the 
error boundary ij\f{n;R) maxje : {R,n,e) is achievable} for fixed rate R. 
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Result 5. For any quantum channel M, the boundary R*{n;e) satisfies 

R*in; e) > Rt^neAn] e), with RAneAn] e) = lA-M) + {e) + O > (15) 

where the coherent information, IfiM), and its variance, Vf{M), are defined in Theorem 6. 

Note that the bound (15) is tight up to the second order asymptotically for the qubit dephasing 
channel (Result 1). However, the bound does not tightly characterize the achievable region of 
general channels, although we have reasons to conjecture that it does for degradable channels. In 
fact, this bound is a direct consequence of an inner bound due to Morgan and Winter [24] together 
with a second order expansion of smooth entropies in [40]. 


III. NOTATION, INFORMATION MEASURES, AND CODES 


In this paper log denotes the binary logarithm. To express the second order expansion of the 
non-asymptotic quantities we need the cumulative standard Gaussian distribution function 


4>(x) 



dy. 


(16) 


We denote finite-dimensional Hilbert spaces by capital letters. In particular, we use A and B to 
model the channel input and output space, whereas M and the isomorphic spaces M' and M” are 
used to model the quantum systems containing the maximally entangled state to be transmitted. 
We also use A” to denote the n-fold tensor product of A for any n G N. The dimension of A is 
denoted by |A|, we use [A] to denote the set {1, 2,..., |A|}, and fix a standard orthonormal basis 

We use £(A) to denote the set of linear operators on A, 'P(A) to denote the set of positive 
semi-definite operators on A, and 5(A) := {p G 'P(A) : tr(/9) = 1} to denote quantum states with 
unit trace on A. A quantum state is called pure if it has rank one. We write /? <C if the support of 
p is contained in the support of a. For general positive operators p,(J ^ V{A), we dehne Uhlmann’s 
fidelity [43] as 

■■= {\\^/pV^\\y , (17) 


where ||A||i := tr [AXX^ is the trace norm. If one of the states is pure, this expression simplihes 
to F{\fi){fi\,a) = (V’kIV’)- 

We often use subscripts to clarify which Hilbert spaces an operator acts on. Let A' be isomorphic 
to A. Throughout the manuscript we denote the maximally entangled state on AA' by 4>aa' = 
\4>){4’\aa' \4>)aa' ~ |A|“^/^ Z]xe[yi] A)A ® A)A'- ^ general state pA G 5(A), its canonical 

purification is {tp^AA' ~ \A\A~PA®'^A’ \A}aa'- clearly have tiA’iA'AA') ~ where tr^' denotes 
the partial trace over A'. 

Quantum channels are completely positive and trace preserving maps between operators and 
denoted by calligraphic letters. In particular, we investigate channels Ma-^b that map 'P(A) 
to V{B). The Choi state of a quantum channel Na^b is defined using the corresponding non- 
calligraphic letter as Nab = (TA ® Ma'^b){\-^Aaa')- 
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A. Information Measures 


Our asymptotic results are stated in terms of the following quantities. For p G Sil-i) and a G 
V{'H) with p a, Umegaki’s relative entropy [44] and the quantum relative entropy variance [21, 
40] are given by 


D{p\\a) :=tr[p {log p-log a)] and V{p\\a) 


tr 


p (log p - log cr 


D{ph)? 


(18) 


respectively. The conditional entropy and the conditional entropy variance [40] of a state pab £ 
S{AB) are given as 


H{A\B)p := -DipABWl-A® Pb) and V{A\B)p := V{pab\\1-A ® Pb) , (19) 

respectively. Related to this we define the coherent information of pab as I{A)B)p := —H{A\B)p 
and its corresponding variance V{A)B)p := V{A\B)p. 

For our non-asymptotic results, we require the following quantities. For p G S{l-L) and a G 'P{T-L) 
the hypothesis testing relative entropy [45] is defined as 


l|/(/5||cr) ;=-log/3i_e(p||cr) with f3i-e{p\\a):= ^mjn tr[A(T]. 

tr[Ap]>l—£ 


( 20 ) 


The hypothesis testing Rains relative entropy of a quantum channel Ma^b is defined as (fol¬ 
lowing the generalized divergence framework discussed in [42]), 

I%{Ma^b)-= sup I%{A : B)j^ ) with I%{A:B)p:= inf Dj^{pAB\WAB), 

Pa(£S{A) aAB&PPT'iA-.B) 

( 21 ) 


where PPT'(A:il) is the Rains set [1, 29], a superset of the set of positive partial transpose (PPT) 
states. It is defined as 


PPT'(A:R) := {tab G r{AB) | \\Tb{tab)\\, < l} , 


( 22 ) 


where Te denotes the partial transpose map on B. In particular, we have the following inequal¬ 
ity [28, Lm. 2]. For every uab G PPT'(A : B), we have 


{4’Wab\(P)ab < 


1 

W\ 


(23) 


for all maximally entangled states \4’)ab of local dimension \M\. Finally, a quantum channel A/a_s._b 
is called PPT preserving if a PPT state input necessarily results in a PPT state output. It turns 
out that PPT-preserving channels output PPT states for any input, since they have PPT Choi 
states [29] (see the discussion after Eq. 4.13). Channels with PPT Choi states were also called 
PPT-binding in [19]. 


B. Codes for Entanglement Transmission Assisted by Classical Post-Processing 

We have defined unassisted entanglement-transmission codes in Section II and in Figure 1. Let 
us reintroduce them in the context of codes assisted by classical post-processing. 

For this, we again consider any quantum channel M = Na-aB and its n-fold extension 
that maps states on A” to states on B'^. An entanglement transmission code assisted by classical 
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FIG. 5: Coding Scheme for entanglement transmission over n uses of a channel JVa^b with classical post¬ 
processing. The encoder S = Sm'^a^q encodes M' into the channel input systems and a local memory Q. 
Later, the decoder T) = recovers the maximally entangled state from the channel output systems 

and the memory Q using classical communication and local operations. The performance of the code is 
measured using the fidelity F{4>mm', Pmm")- 


post-processing for AA®” is given by a triplet {\M\,S,'D}, as depicted in Figure 5. Here, |M| is 
the local dimension of a maximally entangled state \(P)mm' transmitted over AA®”. 

The encoder Sm'^A’^q is a completely positive trace-preserving map that prepares the channel 
inputs Ai,A 2 , ... An and a local memory system, which we denote by Q. The decoder Vqb^^m" 
is a completely positive trace-preserving map that is restricted to local operations and classical 
communication with regards to the bipartition Q : B'^ and outputs M" on the receiver’s side. 

Finally, we note that unassisted codes are recovered if we choose Q to be trivial. Hence, 
unassisted codes are contained in the set of assisted codes. 


IV. PROOFS: GENERAL BOUNDS 

It is convenient to first discuss the general bounds discussed in Results 4 and 5. 


A. General Outer Bounds on the Achievable Region 

In this section we derive the general outer bounds from Result 4, precisely stated as Corollary 2 
and Theorem 3 below. Our results are inspired by the strong converse results for generalized 
dephasing channels from [42] and the metaconverse for classical channel coding [27]. 

We first formulate a general metaconverse bounding possible rates R given a tolerated error e 
for single uses of a fixed channel A/". This bound has the useful property that channel symmetries 
can be used to simplify its form. Nevertheless, when applied to n instances of A/", the bound is not 
efficiently computable. Loosening the bound produces a more computationally tractable convex 
optimization, specihcally a semidehnite program. 

Lemma 1. Let Ma^b be a quantum channel. Then, for any fixed e G (0,1), the achievable region 
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with cpp-assistance satisfies,"^ 

i?^PP(l;£)</|j(AA). (24) 

Proof. First, observe that the encoding operation Sm'^aq can be chosen to be an isometry without 
loss of generality, because we may include any extension systems needed for the Stinespring dilation 
into Q. Then we may express the entanglement fidelity as follows 

F = irlcpMM'T^BQ^M' ° o £M'^AQ{(t>MM')] (25) 

= ^^[£m^AQ ® T^BQ^M'i^MM'WA^BiSM'^AQ ® ^M^AqC'/’MM'))] • (26) 

Since £ is an isometry, the state Paaqq — £m'^aq ® £m^aq{4>mm') is pure, and therefore there 
exists an isometry Wai^aqq such that \p) aAqq = ^A'^AQQ \'4’^)aa'- Thus, 

F = ® '£'^BQ^M'i^MM'))WAi^AQQ-^A^B{i^AA')\ ' 

Now consider the entanglement fidelity of any PPT' state (Ta'b instead of AfA-^Bi'f’AA')- (^^) 
we have 

tT[(l)MM'{£lj^AQ ® '^bq^m'{Wa,^aqq(^a'bW[,^^q^))] < — , (28) 

as the operations on aA'B are all PPT-preserving. We may write this bound in terms of the 
hypothesis-testing relative entropy, because 

Aa'b '■= w\/^aqq^^ai^aQ ® '£'^bq^m'('^aim'))Wa/^aqq ( 29 ) 

is a feasible test to discriminate between A/’a_s.b(V’^/^)) and (Ta'b- That is, Aa'b satisfies 0 < 
Aa'b < Ia'B and tr[AA'BA/’A_5>s(V'A'A))] ^ 1 “ former since V is completely-positive and 

trace-preserving and £ and W are isometries, the latter by assumption that F > 1 — e. From (28) 
we then obtain 

i?=PP(l; £) < D%{Na^b{^%a)) II ^a>b) • (30) 

Since the bound holds for all PPT' aA'B: we may take the infimum over this set to obtain 

(31) 

This bound depends on the precise channel input pA £ ‘5(A) used by the code, but we can remove 
the dependence by taking the supremum over all possible inputs. The result is (24). □ 

Applied to the channel AA®” we immediately get for any fixed e G (0,1), 

R^P'^{n;e) . (32) 

This bound forms a counterpart to Lemma 4, but suffers from the same weakness. It is generally 
hard to evaluate this bound even for moderately large n. However, we may relax the bound 
from Lemma 1 to a convex optimization by restricting the form of the possible states aAB iu the 
definition of the hypothesis testing Rains relative entropy I^. 


Indeed, the outer bound also holds for coding schemes with (unphysical) positive partial transpose (PPT) assistance 
and this includes in particular classical pre- and post-processing assistance (see, e.g., [20] for a precise definition 
of PPT assisted codes). 
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Corollary 2. Let Na-^b be a quantum ehannel. We define the function 

e) := inf inf sup , (33) 

Pa&S{A) Kab&{paWa) Ma^b&PPT 

with the set T{pA,M,e) := {Aab ■ 0 < Aab ^ Pa® ^b A^^^IAabN ab] > 1 — s}, and the Choi states 
Mab of M.a^b and Nab of Ma^b- Then, for any fixed e G (0,1), the achievable region satisfies 

R-PP(l;e)<-logf{J\f,e). (34) 

Proof. Suppose that cjab = iTA®M.A'^B){'f’AA') some PPT-preserving (PPT-binding) channel 
Ma^b- Writing out the right-hand side of (24) using this relaxation and notation gives 

^''PP(l;e) <-log inf sup inf tv [A'j^^^{Ia ® MA'^B)ifi>AA')] ■ (3^) 

Pa&S{A) Ma^b^PP^ 

tr[A^p]>l—e 

Now we may define Aab = (Pa)^'^^^ab(Pa)^^^ 

^'^•’^(1; e) < — log inf sup inf tr[AABM^^] . (36) 

PAeS(A) Ma^bsppt AAser(pA,A/',£) 

Finally for fixed channel input pA, we can reverse the order of the inner optimizations in (36) 
by von Neumann’s minimax theorem, since the objective function is linear and the sets are both 
convex and compact. This concludes the proof. □ 

Furthermore, we show in Appendix A that f(N, e) can be expressed as a semidefinite optimiza¬ 
tion program that satisfies strong duality. 

Group Covariant Channels: In the following we show that symmetries of the channel can 
further simplify the outer bounds. First let us state precisely what we mean by symmetries. 
Suppose G is a group represented by unitary operators Ug on A and Vg on B. A quantum channel 
Ma-aB is covariant with respect to G when 


VgM{-)V^g = N{Ug • C/t), V5 G G. (37) 

Alternatively we can also write this as an invariance of the channel 

M{-) = V^gM{Ug • Ul)Vg, V5 G G . (38) 

Now the main workhorse to simplify our outer bounds for channels with symmetries is [42, Prop. 2], 
which states that we may restrict the optimization in Lemma 1 to covariant input states. Due to 
the form of the hypothesis testing Rains relative entropy, we may then also choose group covariant 
PPT' states a and test operators A to obtain the tightest bound. Note that the convex optimization 
outer bound in Corollary 2 inherits these symmetry simplihcations. 

For general tensor product channels, which are invariant to permutation of the inputs and 
outputs, this allows us to restrict attention to pure states that are permutation invariant. Moreover, 
if the channel is covariant, that is, covariant with respect to the full unitary group, then the channel 
input state can be chosen to be maximally mixed. 

Now let Na-aB be a covariant quantum channel and (f)AA' a maximally entangled state. Then, 
we bound 

R^PP(n;£)< min , (39) 

(TABePPP'(A-.B) n ^ 
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where we voluntarily restricted the minimization to product states in PPT'(74:i?). Moreover, 
since these states have tensor product structure, the outer bound can be expanded using [21, 40] 

= D{p\\a) + + o . (40) 

n V n \ n J 

This leads to the following theorem. 

Theorem 3. Let M = Ma^b be a quantum channel. We define its Rains information as 

■= min D{J\fA'^Bi4>AA')\\crAB) ■ (41) 

<7ab&PPT'{A-.B) 

and let 11 C PPT'(A : B) be the set of states that achieve the minimum. The variance of the 
channel Rains information is 

{ maxV{J\fA'^B{4>AA’)\WAB) for e<4 
-Asen _ ^42) 

m.mV{AfA'^B{4>AA’)\\crAB) for e>2 

CTAB Gil 

If J\f is covariant, then for any fixed e G (0,1), the achievable region with cpp-assistance satisfies 

IfiPP{n]e) < RZterin;£), with RZZin;£) = • (43) 

Since we are here interested in outer bounds, we are also free to chose a potentially sub-optimal 
(^AB £ PPT'(A : B) to further relax this bound. As we will see in Section V for the qubit dephasing 
channel and the erasure channel with classical post-processing assistance the bound from Theorem 3 
is tight up to the second order asymptotically. 


B. General Inner Bounds on the Achievable Region 


In this section we derive the general inner bound from Result 5, stated as Theorem 6 below. 
We use the decoupling approach [13, 14, 18], and in particular a one-shot bound by Morgan and 
Winter [24] which is a tighter version of previous bounds [4, 6]. 

To reproduce their result we need the following additional notation. Sub-normalized quantum 
states are collected in the set 5,(A) := {p G P(A) : tr(/3) < 1}. The purified distance [39] s- 
ball around p G S{A) is then defined as B’^{p) := {p G S,(fH) \ F{p,p) > (1 — e)^}- Finally, for 
Pab £ S{AB) and e > 0 the smooth conditional min-entropy [30, 39] is defined as 


ue 

'^min 


{A\B)p-.= sup sup sup ■ 
Pab&I3^{pab) CTseiS(B) 


A G 


PAB <2 ^ ■ 1a0 O'B 


}■ 


(44) 


Let us now restate Morgan and Winter’s result expressed in terms of the non-asymptotic achiev¬ 
able region as introduced in Section II. 


Lemma 4. [24, Prop. 20] Let Ma-aB be a quantum channel with complementary channel 
Then {R, l,e) is achievable if, for any p G (0,e] and any state pA G So{A), we have 

i?<i^;f""(A|i?).-41ogI, (45) 


where ujae = (2)4 ® A/'|/^£;)(V’( 4 A')- 
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Note that Morgan and Winter use the purified distance as their figure of merit whereas we use 
the fidelity criterion (1). This accounts for the square root in the smoothing parameter of the 
conditional min-entropy. Also Morgan and Winter state their result for the special case n = 1, but 
this can be generalized to arbitrary n G N if we simply consider as a single channel. This 

leads immediately to the following inner bound on the achievable region: 

Corollary 5. Using the notation of Lemma 4 with oja^e^ = ® we have 

R{n;e)> sup sup - - Alog^ - l) . (46) 

vme)PAn&S{A")n'- V 

The problem with this bound is that it is generally hard to evaluate, even for moderately large 
values of re. Hence we are interested to further simplify the expression on the right-hand side in 
this regime. To do so, we choose rj = \l\fn and use input states of the form p®”. This yields the 
following relaxation, which holds if re > 

R{n;e)> sup - 21ogre - l) . (47) 

pa&s{A) n 

Here we introduced Sn = y/e—)= and ujae as in Lemma 4. Using standard second order expansion 

■y ^ 

methods [40], we can give an asymptotic expansion of i?*jj^gj,(re; e) in (47) as follows. 

Theorem 6 . Let N = Ma^b be a quantum channel. We define its coherent information as 
IciAf) := max I{A)B)^, with lvab = (d^A ^ AfA'^B){fi’AA') 

pa&S{A) 

and let H C 5o(A) be the set of states that achieve the maximum. Define 


VAN) := 


rminl/(.4)B). 

I max V(A)B)i, 

\pa&ti 


for 

e < 2 

for 

^>5 


Then, for any fixed e G (0,1), the achievable region satisfies 


(49) 


R{n]e) > Anner(n;e), with Rir,ner{n] e) = Ic{M) + \j '' ^(e) O . ( 50 ) 

V re \ n J 

Proof. We analyze the expression in (47) using the following asymptotic expansion of the smooth 
conditional min-entropy [40], 

= H(A\B), + $-1 (j2) + o (!^) . (51) 

n ^ \ n \ n J 

This yields that for any pA G S{A)^ we have 

R{n-,e)>H{A\E)^ + J^^^^^-\e) + o(^^) , (52) 

\ n \ n J 

and then by duality of the conditional entropy we find H{A\E)i^ = I{A)B)^. Furthermore, it is 
easy to verify that V{A)E)i^ = V{A)B)i,j (see, e.g., [17]). We conclude the proof by choosing an 
optimal state pA^Li depending on the sign of <^“^( 6 ). □ 


As we will see in Section V for the qubit dephasing channel the bound in Theorem 6 agrees 
with the outer bound stated in Result 4 up to the second order asymptotically. 
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V. PROOFS: EXAMPLES 

In this section we derive our results concerning the qubit dephasing channel, the erasure channel, 
and the qubit depolarizing channel as announced in Results 1, 2, and 3. 


A. Covariant Generalized Dephasing Channels 

First we consider covariant generalized dephasing channels, which have the (additional) property 
that has full support on the projector 11 = ( 8 ) \x){x\b in some basis. 

In that case, starting from (39), we can use the data-processing inequality for a map S{-) = 
n • n + (1 — n) • (1 — Ft) to write 

R^PP(n;e)< min -D^j,{{ZA'^B{cl>AA'))^^h%) (53) 

< min (54) 

ugeSiB) n D / 

< izi|,((ZA^^s((/)AA0)®1|(lA®^A'^s(<^A'))®") • (55) 

To show the second inequality we apply T®” to both states and employ the data-processing in¬ 
equality for the hypothesis testing relative entropy. Note that the map £ keeps Z{(j)AA') invariant 
and maps 1 ^ ( 8 ) to a normalized state aAB that is classically correlated and thus in particular 
in PPT'(A : B). The bound in (55) has the form of a conditional entropy or coherent information, 
and it can be expanded again using (40) to find 

R"PP(n; e) < + o . ( 55 ) 

So for covariant generalized dephasing channels this now agrees with the inner bound from The¬ 
orem 6 up to the second order asymptotically (an example being the qubit dephasing channel as 
discussed in Section VB). Given the recent results [42] the outer bound (56) might also hold for 
generalized dephasing channels (without assuming covariance) and then agree with the inner bound 
from Theorem 6 up to the second order asymptotically. 


B. The Qubit Dephasing Channel 

The qubit dephasing channel is defined as 

Z^ : {l-j)p + -fZpZ , (57) 

where 7 G [0,1] is a parameter and Z is the Pauli Z operator. This channel is covariant since it is 
a qubit Pauli channel. Now to determine the second order asymptotic performance it is sufficient 
to specialize the outer bound from (56) and to apply Theorem 6 for the inner bound. It is then 
easily seen that 


= with /i(7) = 1-ylogy-(1-7)log(l-7) 

V{R)B)z.y{p) =v{'y) with u(7) = 7(log7 -F /i(7))^-F (1 - 7)(log(l - 7 ) + ^( 7 ))^ 


(58) 

(59) 
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and hence we deduce 

R{n-,e) = l-h{j) + \[^<^-\e) + o(^-°^) . (60) 

V n \ n J 

However, we can refine (60) and determine the third order asymptotic performance. We do this 
by directly obtaining the finite block length behavior of the qubit dephasing channel from that of 
the classical binary symmetric channel (BSC). First, consider the converse, particularly that of (39), 
applied to the channel Z®'^. Using the Bell states 4)^^ = (I)aa' and 4’'^^ = {1a®Za')4>AA'{^a®Za')^ 
we immediately find 


UJAB ■= Ma'^b{4aA') = (1 - l)(t>AB + ^^AB ' ( 61 ) 

Now, in (39) we are free to pick any PPT' state to obtain a bound. Pick aAB = ^ab)^ 

which gives^ 

R{n-,E) < < \iy„(^%\\cT%) . (62) 

To connect to the finite block length bounds of the BSC, consider measuring both A and B in 
the Pauli x basis, and let X and Y be the output random variables for A and B, respectively. 
For the state ujab, this results in the distribution Pxy in which Px is uniformly-distributed and 
P[Y = X] = 1 — 7 . For aAB, the distribution is of product form PxQy with Qy also uniform. 
Moreover, the original quantum states can be reconstructed from the classical random variables 
X and Y by the map which outputs 4~ab '''^hen X = Y and 4>~ab otherwise. Therefore, the bound 
becomes 

R{n-e) < x , (63) 

which is precisely the bound obtained by Polyanskiy et al. for the BSC [27, Theorem 26], which 
is equivalent to the classical sphere-packing bound [15, Eq. 5.8.19]. This establishes one inequality 
(upper bound) in ( 6 ). 

For the achievability, we may directly employ linear codes for the classical BSC to the qubit 
dephasing channel. Specifically, any linear {R, re, e} code for the BSC (which recovers the input 
with probability at least 1 — e, averaged over a uniform choice of inputs), can be converted into 
an {R, re, e} Calderbank-Shor-Steane (CSS) code for entanglement transmission over the dephasing 
channel. This is possible since, for a linear code, the action of the channel is a mapping among the 
orthogonal Bell states, which is essentially a classical action. 

To formalize the connection, begin with the description of the classical linear code by its (log re — 
log|M|) X re parity check matrix H. Each row rj G {0,1}"' defines a parity function and the 
codewords Ck of the code must satisfy Ck ■ rj =0 for all j. The associated CSS code can be 
defined as the simultaneous -|-1 eigenspace of the “stabilizer” operators X^'R where Z 

■ ■ ■ Z Crucially, the action of the channel is to apply an operator of the form with 

u G {0,1}", according to the distribution Pu. At the output, the receiver can simultaneously 
determine the eigenvalues of all the of the stabilizer operators. This information is precisely 
equivalent to determining the value of the parity checks of the classical linear code, called the 


® The choice of aAB is equivalent to using the convex relaxation of the bound, Corollary 2, and choosing M = 
in (33). 

® Generically, a CSS code has stabilizers of both X-type, as here, and of Z-type, i.e. composed of products of Pauli 
Z operators. 
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syndrome s. Given the syndrome, the decoder of the classical code determines a guess as to the 
input codeword, which is equivalent to a guess u'{s) of the actual channel error. 

We may also utilize this algorithm (whatever its precise details) in the quantum case, and 
attempt to correct the error by applying When u'{s) is the true error pattern, the quantum 

state is properly recovered, and the entanglement fidelity is unity. On the other hand, if u'(s) is 
incorrect, then in the worst case the action (^)+“ is a logical operation on the code subspace, 
which results in a state orthogonal to the desired entangled state. Therefore, the error probability 
of the classical code translates directly into the entanglement fidelity of the quantum code. Thus, 
we may apply finite-block length bounds for linear codes, particularly the bound of Poltyrev [25] 
(see also [27, Eq. 65]). This establishes the other inequality (lower bound) in (6). 


C. The Erasure Channel 

The qubit erasure channel is defined as 

Sis : {1 - I3)p +/3\e){e\ , (64) 

where /3 G [0,1] is a parameter and le)(ej is a quantum state orthogonal to p. Using the covariance 
of the channel, we could first obtain a second order asymptotic similar to that of the dephasing 
channel in (60). However, it is not too difficult to directly derive an outer bound and an explicit 
coding scheme leading to an inner bound, which precisely match for all n. 

Let us begin with the outer bound. Again we may relate the finite block length performance 
to a classical coding problem, namely the classical binary erasure channel (BEC). The argument 
for the outer bound proceeds very similarly to the dephasing example. The optimal channel inpnt 
state corresponds to the maximally entangled state and the state prodnced by the channel 

is now 


AB = (1 - P)4>ab + Ptta G [e)(el^ , (65) 

where tta denotes the maximally-mixed state. Measurement of A in the Pauli x basis and B in the 
basis {]+),]—), Je)} produces the distribution Pxy with Px uniform and Y = X with probability 
1 — /3 and Y = e with probability 1-/3. The original state can be reconstructed using the map 
which sends {X, Y) to when X = Y, when X / T / e, and to tt^ ® \^){^\b when Y = e 
otherwise. As before, we make a specific choice of PPT' state in (39), but this time not a product 
state. Instead, consider the classical distribution P^"' x Qy^ given in [27, Eq. 168]. The Qyn 
distribution has the property that any two y” with the same number of erasure symbols e have the 
same probability, i.e. there is no dependence on the number of Os versus Is. The aforementioned 
map takes the distribution to a quantum state which is diagonal in the standard bases {jO) , jl)} 
for A and {jO) , jl) , je)} for B, and is therefore a PPT state. This can be seen as follows. Consider 
a fixed position j in a given a pair (x"^, y"). If yj = e, the state of the jth pair of systems AB is 
manifestly diagonal in the standard basis. On the other hand, if yj ^ e, then the state is mapped 
to either 0]^^ or depending on the value of xj. But the sequence in which yj takes the other 
value has identical probability, meaning the two Bell states occur with equal probability, making 
the AB state diagonal. Since we may map and aA" b^ to the associated classical distributions 
and back, the following converse holds for the qubit erasure channel, 

fppv^n-e) < >< Qy-) • 


( 66 ) 
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By design in the choice of ctabJ this is precisely the bound for the BEC reported by Polyanskiy 
et al. [27, Thru. 38], as discussed in more detail by Polyanskiy [26]. 

Next, we construct an explicit coding scheme, involving classical post-processing including com¬ 
munication from the receiver to the sender, which matches the outer bound exactly. (See Figure 5 
for schematic description of this code.) The strategy of the coding scheme is to generate maximally 
entangled qubit states using the quantum channels and then use the successfully transmitted (i.e. 
not erased) maximally entangled qubit states to distill a an entangled state of local dimension |M|, 
as required. (Note that the number \M\ is fixed at the outset of the code, i.e. the entanglement 
transmission scheme must deliver a maximally entangled state with local dimension |M|, possibly 
at the expense of low fidelity, rather than outputting a variable number of certifiably high fidelity 
entangled pairs.) 

First, the encoder prepares n maximally entangled qubit states IV’) and sends one half of each 
over the channel. The other halves, together with the untouched system M', are stored in the 
memory register Q. The decoder now works as follows. The receiver determines which qubits have 
not been erased and informs the sender of their locations. Let L be the random variable indicating 
the total number of erasures and note that L follows a binomial distribution with parameters n 
and /?. Let us also fix fc = [ log |M|] and consider the following two cases: 

1. If L = /<n — A; the decoder can extract a maximally entangled state with unit fidelity. To 
do so, it selects k perfectly transmitted entangled qubits at the sender and receiver. Let us 
assume (without loss of generality) that these are in a state K*)- 

The receiver then prepares a maximally entangled state of local dimension \M\ by measuring 
the k qubits with the projective measure 

[2‘]^l‘S| = l"l[ ■ (67) 

[ i&S ) 

The outcome, a subset S of cardinality \M\, is transmitted to the sender so that both sender 
and receiver now share a maximally entangled state on the subspace determined by S. 


2. On the other hand, \i L = I > n — k sender and receiver simply select the successfully 
transmitted qubits and embed them in a space of local dimension \M\. The fidelity with the 
target state \(f)) = —?= N) is given by 


VWi 




|M| 




~)n—l 


\M\ 


( 68 ) 


To complete the decoding operation, the sender and receiver perform quantum teleportation to 
teleport M' to the receiver, using the maximally entangled state prepared above as a resource. The 
fidelity of the state prepared above with the target state 4 > mm ' is then just the expected fidelity 
over L, which evaluates to 

7^ = "f (>'(!-«”-'+ t (>'('- 

1=0 ^ ^ l=n-k+l ^ ' 


This corresponds to using Corollary 2 with M the channel which ignores its input and prepares (jb-^ at the output. 
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This is exactly the expression reported in the aforementioned outer bound in [27, Thm. 38], mean¬ 
ing the inner bound coincides with the outer bound when we allow classical post-processing and 
communication from the receiver to the sender. 


D. The Qubit Depolarizing Channel 

The qubit depolarizing channel is defined as 

ct 

V^-.p^{l-a)p+-{XpX + YpY + ZpZ) , (71) 

where a G [0,1] is a parameter and X, T, Z are the Pauli operators. This channel is covariant since 
it is a qubit Pauli channel. Using the Bell states (Pab ~ ® Za')(Paa'{^a ® Za'), 

^\b = 0-A ® Xa’)(I)aa'{'^a ® Xa'), and = {1a ® YA')(pAA'ilA ® Pa/), we immediately find 

Q/ 

ujab ■= (Pa ® VoP){(pAA') = (1 - a)0As + 3 ('/'as + '^ab + 'P’ab) ■ ('^2) 

Now choosing aAB = \(P\b P \ ^^~ab P ''/’as P ^ab) ™ (^^) the outer bound 

< lofeMSIkSfl) ■ (73) 

As in the case of the qubit dephasing channel, we can convert the hypothesis test between ojab 
and aAB into a test between classical distributions, in fact precisely those distributions which were 
used in the dephasing example. This follows by considering the map which generates (p\^ when 
X = Y and otherwise randomly generates one of the other Bell states when X ^ Y. Therefore, 
we obtain the same outer bound for the qubit depolarization channel as for the qubit dephasing 
channel. This raises the question of whether cpp assistance (or more generally PPT assistance, 
cf. Footnote 4) can turn the qubit depolarizing channel into the qubit dephasing channel. 


VI. CONCLUSION 

We gave inner (achievability) and outer (converse) bounds on the boundary of the achievable 
region for quantum communication with finite resources. We showed that these bounds agree for 
the qubit dephasing channel and qubit erasure channel with classical post-processing assistance 
up to the second order asymptotically. Moreover, we even gave a third order characterization for 
these specific examples by employing finite block length bounds of the binary symmetric channel 
and the binary erasure channel, respectively. 

However, many questions remain open. For example, we would like to understand if the inner 
bound in Result 5 characterizes the achievable region for all degradable channels [10] (cf. the open 
questions in [42] ). Also it would also be interesting to explore higher order refinements for channels 
with zero quantum capacity (e.g., for the erasure channel with /3 > 1/2 and no assistance). This 
might lead to a better understanding of super activation of the quantum capacity [36]. 

Finally, we would like to note that the recent results about finite resource entanglement assisted 
classical communication [8] can immediately be transformed to entanglement assisted quantum 
communication (and this then also gives outer bounds on the achievable rate region for unassisted 
codes). This is accomplished by using the equivalence results in [20, App. B] which make use of 
quantum teleportation and superdense coding. In particular, one finds that for covariant channels 
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unassisted achievable region 
unassisted capacity 

entanglement-assisted achievable region 
entanglement-assisted capacity 


FIG. 6: Second order approximation of the achievable region of a qubit dephasing channel with e = 1% and 
7 = 0.1 ; the achievable region is enlarged in the presence of entanglement [8]. 


J\f (which includes the qubit dephasing channel and the erasure channel) the boundary of the 
entanglement assisted achievable region i?^(n;e) satisfies 


R^{n] e) 


2 


4n 



(74) 


with the mutual information of the channel, I{M), and its variance, V^{J\f), as defined in [8]. As 
an example, we mention again the qubit dephasing channel for which 


i?^(n; s) 


1 - 2h(7)+ 



^-^e) + 0 



(75) 


where h{'y) denotes the binary entropy and v{'y) the corresponding variance as defined in Result 1. 
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Appendix A: Semidefinite Optimization 

In this section we describe how to formulate the outer bound from Corollary 2 as a semidefinite 
optimization program satisfying strong duality. 

A semidefinite program (SDP) is simply an optimization of a linear function of a matrix or 
operator over a feasible set of inputs defined by positive semidefinite constraints. We give only 
the bare essentials here, for more detail see [5, 46]. The maximization form of an SDP is defined 
by a Hermiticity-preserving superoperator £a^b taking C{A) to T(R), a constraint operator C G 
C{B), and an operator K G Z{A) which defines the objective function. The SDP is the following 
optimization, which we will also refer to as the primal form, 

a = supremum tr[AA] 
subject to £{X) < C 
X >0. 


(Al) 
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When the feasible set is empty, i.e. no X satisfy the constraints, we set a = —oo. The dual form 
arises as the optimal upper bound to the primal form, and takes the form 

(5 = infimum tr[C'y] 

subject to S*{Y) > K (A2) 

y > 0. 


Again, when the set of feasible Y is empty, /3 = oo. Weak duality is the statement that a < /3, 
that indeed the dual form gives upper bounds to the primal (or that the primal lower bounds the 
dual). Strong duality is the statement that the optimal upper bound equals the value of the primal 
problem, a = /?. This state of affairs often holds in problems of interest, and can be established by 
either of the following Slater conditions. In the first, called strict primal feasibility, strong duality 
holds if /3 is finite and there exists an A > 0 such that £{X) < C. Contrariwise, under strict dual 
feasibility strong duality holds when a is finite and there exists a T > 0 such that £*{Y) > K. For 
strongly dual SDPs we also have the so-called complementary slackness conditions £*{Y)X = KX 
and £{X)Y = CY that relate the primal and dual optimizers. 

Theorem 7. With the notation from Corollary 2, the outer bound f{N, e) can he written as 


f{J\f,s) = infimum tr[(^yi] 

subject to Ca ® > Aab + ©as 

^AB £ £{pa-,J^ )£)j Pa £ ‘5(A) 
Ca, Qab > 0, 

or, equivalently, 

f{M,e) = supremum — e) — v 

subject to pNab < Mab + Rab 
tlCB[RAB] < nlA 
Mab G ppt; p, v, Rab > 0. 


(A3) 


(A4) 


Proof. The proof is straightforward: we simply use the dual of the inner optimization in (33) to 
obtain the minimization problem (A3). Then we use Slater’s condition to show that strong duality 
holds and obtain (A4). 

Consider the function 


MOab) ■= sup tT[OABMAB] , (AS) 

yV( a->sSPPT 

and observe that /o is a semidefinite program. In particular, it is a primal problem as we have 
defined it, with X = Mab, K = Oab, C = (0, 1a), and £{X) = {—X'^^,iTB[X]). Choosing for the 
dual variables Y = (Tabj^a), the dual of /o is 

MOab) ■= infimum tr[^A] 

subject to Ca ® Is > Oab + ^'ab 
£ab,^a > 0 . 

Combining this with the outer optimization in (33) gives the minimization program (A3). The 
equality statement is precisely strong duality of the primal and dual forms of the inner optimization. 
By Slater’s condition, strong duality holds if /o is finite and there exists a strictly feasible set of 
dual variables. Observe that /o(Oas) < |A|, since for the optimal Mas we have /o(Oas) = 
tr[MAsOAs] < tr[MAs] < = |A|. Here we have used the upper bounds Oab < Ias and 
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ii b[Mab\ < 1a- Thus, the first condition is fulfilled. Meanwhile, Tab = Ias and = 3 • are 
a strictly feasible pair. Thus, /o = /o over the domain of interest. 

To construct the maximization program, we simply dualize the minimization program. In 
particular, f{Af,e) is a dual-form semidefinite program in the variable Y = {4>a,A.ab,^ab,Ca) 
with C = (0,0, 0, Ia), K = {1 — e, —1, 0, 0), and 

£*(Y) = {ti[NAB-^AB],- tr[(/>A], (/'a ® Is - ^ab,Ca ® Is - ^ab - T^b) • (A7) 

Choosing primal variables X = {m,n, Rab, Mab) leads to the maximization in (A4). Equality 
again follows from Slater’s condition: / is finite (in particular the bound on /o used above), while 
a feasible choice of dual variables is given by Mab = Rab = ^^\^ab-, s = 1 , and m = 2 \a\\b\ • 
choice of m ensures the first constraint holds strictly, since any Choi operator of a trace-preserving 
map satisfies ||A^as||oo = |A| (largest singular value). □ 

No discussion of strong duality of semidefinite programs is complete until the complementary 
slackness conditions have been formulated. Often, these give considerable insight into the form 


and properties of the optimizing variables. First observe that 

£{X) = {-uIa + tiB[R^AB\^ rnNAB - Mab - Rab, ItbIMab]) ■ (AS) 

Then the conditions are easy to read off from the form of C and K. They are 

ix[(f)A] = 1 (A9) 

tTc[AABNAB] = l-£ (AlO) 

4>aRab = AabRab (All) 

^aMab = {Aab + '^'^ab)Mab ( A 12) 

n(j)A = tvB[RA%](t)A (A13) 

Kb^ab = 0 (A14) 

trB[MAB]^A = '^A (A15) 

vaNabAab = {Mab + Rab)Aab ■ (A16) 
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